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SUMMARY 


Viscous self-induced pressures on 3^-semiver-tex-angle cones were 
measured over the range 3. 7 < free-stream Mach mmiber <5.8 and 0. 5 < 
viscous-interaction parameter < 2.3. The data were found to he in good 
sigreement with results obtained by Talbot on 5^^ cones in the range 
3. 7 < free-stream Mach n-umber <4.1 and 0. 9 < viscous -interaction pa- 
rameter < 3. 5. All these data were correlated reasonably well by the 

viscous-interaction parameter, which is defined as = Mg -y/C/Re^^, 
where and "the Mach number and Reynolds number based on 

ideal Taylor -Mac coll flow conditions and C is the Chapman-Rubesin 
factor. 

A new method for calculating self-induced pressures is presented 
which takes into accouint the interaction between boimdary-layer growth 
and the inviscid-flow field at the outer edge of the boundary layer. 
Pressiires calculated by this method were only 10 to 20 percent higher 
than the measured values. 


INTRODUCTION 

The fluid-dynamic and thermodynamic phenomena associated with flight 
at hypersonic speeds have been the subject of intensive research in re- 
cent years. In this research one problem that has received considerable 
attention is the problem of the "self- induced-pressure" effect, which is 
one aspect of a broader class of phenomena that can be described as 
" viscous- interaction" phenomena. 

The self- induced-pressure effect comes about in the following way. 

At hypersonic speeds the boundary layers which develop on bodies are, 
because of the large temperature differences generated through them, many 
times thicker than those which are produced at low speeds. Since the 
density of the hot gas in a hypersonic boxmdary layer is very low, the 
mass flux within the bouindary layer is small. Thus, the presence of a 
thick layer of hot gas adjacent to the surface of a body results in the 
outward displacement of streamlines in the flow external to the layer; 
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this outward displacement can he regarded as equivalent to an effective 
thickening of the body (see fig. l). It can easily be seen that this 
thickening will result in increases in pressiure in the flow external to 
the boundary layer, and because the pressure in the external flow is 
transmitted essentially without change through the boundary layer, the 
pressures on the surface of the body will likewise be increased. The 
difference between the actvial surface pressure and that calciilated by 
inviscid theory neglecting the boundary layer is called the self- induced 
pressure. 

Although the fundamental mechanism responsible for self-induced pres- 
sures is well understood, the analysis of the effect is rather compli- 
cated. The magnitude of the self-induced pressirre is directly proportion- 
al to the rate of growth of the boundary layer. However, the growth of 
the boimdary layer is determined by the pressvtre, Mach number, and so 
forth in the flow at the outer edge of the layer, and the values of these 
quantities depend on the magnitude of the displacement effect. It is 
seen, therefore, that the problem under consideration is a complex inter- 
action phenomenon in which the boundary- layer history plays an impor- 
tant role. It will also be recognized that the magnitude of the phenomenon 
is more significant for thin bodies such as flat plates and slender cones 
than for thick bodies, since for thin bodies the changes in effective ge- 
ometry due to boundary- layer growth will be proportionately larger. 

The particular problem considered in this report is the effect of 
self- induced pressxore on a slender cone. In the first part of the re- 
port new data are presented for pressures on 3° semi vertex angle cones. 

The two appendixes disciass several of the analyses which have been pro- 
posed for predicting the effect and present a new method of calculating 
self- induced pressures. 

This investigation was carried out at the University of California 
\inder the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 

C Chapman-Rube sin factor in relation {y-/\i-2) - C(t/T2) 

d diameter of pressinre orifice 

similarity parameter, M2^0(, 

K2 similarity parameter, Mj^02 

M Mach number 

p static pressvire 


L 


4937 


C't- 1 back 


HACA TN 4327 


3 


Re 

Reynolds number 

r 

radial distance measured outvard from cone axis 


cone radius 

T 

gas temperatiare 

t 

bluntness of cone tip or thickness of leading edge of plate 


gas velocity 

X 

distance meas\ired from vertex along cone surface 

y 

distance normal to cone surface 

r 

ratio of specific heats 

o 

boundary-layer thickness 

^ 6* 
rH 

boundary-lay^er displacement thickness, ^ - (pu/p2U22J^ 


displacement thickness on cone ^ corrected for transverse curvature 


cone semi vert ex angle 

r V 


-1 { d6*\ 

streamline inclination at outer edge of boimdary layer, tan 1 1 

®2 

% 

A 

mean free path 


absolute viscosity 

V 

kinematic viscosity 

P 

gas density 

0 

Prandtl number 


•z 

viscous-interaction parameter, M -^JCfRe^ 

^ c 


Subscripts : 


aw 

adiabatic wall 

c 

ideal flow conditions along siorface of cone obtained for x -*■ ” 

w 

wall conditions 

A 

property evaluated at distance A from wall 
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1 fi*ee-stream conditions 

2 conditions at outer edge of cone boundaiy layer ^ taken to be func- 

tions of X 

Examples of the notation used for the various Reynolds numbers are: 

Re 2 / inch = ^ 2 /^ 2 ^ ^^t " and Re^^ = u^x/v^ 

1 c 

PREVIOUS EXPERIMENTAL RESULTS 

The only data which have been available up to now for viscous in- 
teracting flow over a cone seem to be from the experiments of Talbot 
and Baldwin (refs. 1 and 2 , respectively). Talbot obtained his data in 
the Low-Density Wind Tunnel of the University of California. His tests 
were carried out on 5^ -semi vertex- angle cones, in the Mach number range 
3.7 < Mj_ < 4.1, and at Reynolds numbers which corresponded to the viscous- 
interaction-parameter range 0.9 < < 3.5. Baldwin obtained data at the 

Guggenheim Aeronautical Laboratory, California Institute of Technology, 
in the 5- by 5-inch hypersonic wind tunnel, leg number 1. His data were 
also for a 5^ cone, at M]^ = 5.8, and 0.1 < < 1.6. 

APPARATUS 
Wind Tunnel 

The experiments for the present work were conducted in the No. 4 
Low Density Wind Tunnel of the University of California. This wind tun- 
nel, which is an open-jet continuous -flow type employing axially sym- 
metric nozzles, has been described in detail in several publications 
(refs. 1, 3, 4, and 5). Two nozzles were used in the tests: the No. 8 

nozzle, which produces flows in the range 3.7 < < 4.1, 900 < Rej^/inch 

< 3,600; and the No. 9 nozzle, which produces flows in the range 
5.5 < Mj^ < 5.8, 4,000 < Rej_/inch < 9,000. Actual values of the flow pa- 
rameters obtained in the tests are listed in table I. 


Models 

All the cones tested were of 3° semivertex angle. Two sets of 
models were used, each set consisting of seven cones. The type A cones 
had base diameters of 0.500 inch, and the type B cones had base diameters 
of 0.750. inch. The longer type B models were designed primarily to in- 
vestigate the influence on the cone surface pressure of the expansion 
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wave generated at the junctvire of the conical surface and the cylindrical 
afterbody. 

Each cone had four 0.010- inch-diameter press\rre orifices spaced cir- 
cumferentially at 90° intervals at a prescribed distance from its vertex. 
The vertex distances are listed in figure 2. Also shown in this figure 
is the collar and set screw arrangement at the base of the cross-stream 
supports of the models which was used to aline them parallel to the axis 
of the wind-timnel nozzle. Figure 3 shows the type A models, together 
with an impact probe, moimted on the rotary probe selector of the tra- 
verse mechanism. 

Subsequent to its use in the pressure measurement tests, model B7 
was fitted with four copper-constantan thermocouples soldered in the cone 
surface at 2, 3, 4, and 5 inches from the vertex. This model provided 
information on the wall temperature of the cones, which was required for 
the boundary- layer calculations. 

Some difficialty was encountered in producing models with sharp tips . 
The method of fabrication which was found most satisfactory was an acid 
etching process. After the cones had been machined, the tip regions were 
etched by a flat tool covered with a thin film of nitric acid. This 
method produced tips with diameters less than 0.001 inch. However, the 
etching was not completely uniform. In the etched regions of the models, 
which extended back from the vertices about 0.4 inch, some local varia- 
tions in cone angle of several degrees were observed. Farther back from 
the vertices all cone angles were found to be 3.03°±0.03°. 


Instrumentation 

Cone surface press\ares were measured with the temperatixre -regulated 
thermistor manometer described in reference 3. This instnmient was used 
because of its small gage volimie and short tabulation, which results in 
a fast time response, and because of its high accuracy in the measurement 
of the difference between nearly eqml pressures. The thermistor is 
mounted on the traverse mechanism of the wind tunnel, and pressures within 
the models are led to it through a valve mechanism on the rotary probe 
selector. The least count (O.l mv) of the potentiometer used to measure 
the bridge imbalance of the thermistor measuring circuit corresponded to 
a pressure increment of about 0.02 micron of mercury. Before each test 
the thermistor was calibrated statically against a precision McLeod gage 
(ref. 6). Analysis of the thermistor calibration data yielded a probable 
error in absolute pressure of about 1 percent. It took about 1 minute 
for the cone-thermistor system to come to equilibrium after a change in 
pressure was imposed on a model by placing it in the flow. 
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Wind-tunnel stagnation pressinres were measured with a mercury ma- 
nometer. Impact pressiires were measured with a butyl phthalate oil ma- 
nometer (ref. 7). Both manometers are equipped with magnifying optics 
which make it possible to locate the menisci to within 0.001 inch. 


Nozzle Csilibratlon 

Mach numbers and static pressures in the test regions of the noz- 
zles were determined by measuring stagnation and impact pressures, assum- 
ing the flow to be isentropic. For each flow condition of the tests an 
axial impact-pressure survey was made to determine the Mach number and 
static -pres sure variations in the regions occupied by the models. 

The impact-pressure s\irveys in the No. 8 nozzle revealed a region 
about 8 inches in axial extent over which the Mach n\miber variation was 
less than 2 percent and the static-press\are variation less than 10 per- 
cent. In the No. 9 nozzle the axial extent of the region over which 
the Mach number and static -pres sure variations were less than these val- 
ues was about 3.5 inches. 

TEST PROCEDURE 

The test models and an impact probe were mounted on the rotary probe 
selector of the traverse mechanism and allowed to outgas under vacuum for 
2 days. At the beginning of a test the flow was established and measiored 
by means of the impact tube, under "balanced- jet" conditions, which re- 
quired that the pressure measured on the wall of the nozzle be equal to 
the static pressure in the chamber surroimding the open Jet. (This bal- 
ance condition did not alter appreciably when the models replaced the 
impact tube in the stream. ) 

After flow was established, the test models were successively ro- 
tated into the stream, and cone surface pressures were measured. Before 
a measurement was made, sufficient time (about 10 to 15 min.) was allowed 
to permit the model to come to thermal equilibrium; during this period 
the cone pressure decreased slowly because the cooling of the model 
caused a thinning of the boundary layer. 

In the tests reported, all models were positioned so that their 
Y0i.-tices were located at the same point in the flow. Correction was 
made for the axial gradients. Other tests were made with models posi- 
tioned so that their respective pressinre orifices were at the same axial 
location in the stream. The surface pressiares obtained by this second 
method agreed quite well with those obtained by the first method. In 
"the first method the correction for axial gradients in the stream was 
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accomplished sln^jly by using the local Mach numbers obtained from the 
impact-pressure siarveys to determine values of the inviscid surface pres- 
sure Pg corresponding to the actual location of the pressiare orifices. 
These values were used with the measured cone pressiores P 2 to obtain 
the ratio P2/P(,' 

It was fo\md that the above correction procedure gave consistent re- 
s\ilts provided the static pressiire in the flow at the point where orifices 
were located did not differ by more than 10 percent from that at the ver- 
tex of the model. Consequently, in the No. 9 nozzle, where the axial ex- 
tent of usable flow was about 3.5 inches, only models Bl, B2 , and B3 were 
used for the final tests. No difficulties arose on this score in the 
No. 8 nozzle. However, it was observed that the expansion-wave reflec- 
tion of the bow shock (the reflection occurring in the region of the 
strong density gradient where the isentropic core merged with the nozzle 
boundary layer) affected the pressures measijred on models B6 and B7. For 
this reason data are not reported at M = 4 for these models. 

The temperatiire measxnrements made with model B7 showed the surface 
of the cone to be isothermal. The apparent recovery factor (based on 
inviscid-flow conditions Me and T,, behind the conical shock) was 
about 0.89; the increase over the theoretical value of 0.85 for the lami- 
nar boundary layer is due to heat conduction through the support from the 
model mo\mting, which was at essentially stagnation temperature. 


RESULTS 

Results of the tests are given in terms of the induced-pressxire in- 
crement (p 2 /Pc) “ 1- inviscid-flow pressures for the 3° cones were 

calculated from Van Dyke's second-order theory (ref. 8), since interpo- 
lation in the Kopal tables (ref. 9) between 0° and 5° was not suffi- 
ciently accurate.. The data at different flow conditions are shown plot- 
ted against x in figure 4, and in figinre 5 against the viscous- 

interaction parameter = M^ V C/Re^^ , where M(, and Re^^ are the 

Mach n\jmber and Reynolds nxomber based on ideal Taylor-Maccoll flow 
conditions . 

Numerical values for the induced-pressure increment are in the range 
0.06 to 0.30. The acttial measured cone pressures varied between about 
60 and 170 microns of mercury. An error of 1 percent in the measured 
pressijre is equivalent to an error in the induced-pressure increment of 
5 percent or more for most of the data, assuming the values of p^ to 
be exact, from this it is estimated that the over-all probable errors 

in (P 2 /Pc) - 1 between 5 and 15 percent. Free-stream Mach numbers are 
acciarate to about 1 percent; free-stream Reynolds numbers are accurate to 
about 5 percent. 
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In addition to the results of the present tests ^ two sets of data 
taken from reference 1 are shown in figure 6. In figiu*e 5 all the ref- 
erence 1 data for 5° models A1 through A7 are plotted. The data from 
reference 2^ which are also for 5° cones ^ are represented hy a single 
line in figure 5. 


DISCUSSION 
Experimental Results 

One conclusion which can be drawn from an examination of figure 4 
is that the effect of the shoulder expansion did not extend far enough 
upstream in the cone boundary layer to influence the cone surface pres- 
sures. This is evidenced by the fact that the data obtained with the B 
models agree ^ within experimental scatter, with those obtained with the 
A models. It was not possible to determine the extent of the region on 
the cone which is influenced by the shoiilder expansion, because the re- 
flection of the bow shock wave back onto the models obscured the effect. 
Models B6 and B7, which were designed to measure the upstream influence, 
were those most affected by the reflected wave. In figure 4 it can be 
noted that the pressure orifice on model B5 was also within the zone in- 
fluenced by the reflection. 

The scatter in the data is probably due to a combination of experi- 
mental error, imperfections in the orifices, and inaccuracies in the cone 
angle in the tip regions of the models. One may note that the reproduci- 
bility of the data was quite good, as evidenced by the comparisons shown 
in figures 4(e) and 4(f) between different sets of measurements made with 
both the A and the B models. 

It can be seen from figure 5 that the parameter provides a 

fairly good correlation for all the data obtained in the University of 
California Low Density Wind Tunnel. The Mach 4 data for the 3^ and 5^ 
cones agree quite well; the Mach 5.6 data are slightly lower. 

It can also be seen from figure 5 that the induced-pressure incre- 
ments of reference 2 are higher, by a factor of about 2, than those ob- 
tained here. One suggestion which has been advanced is that the differ- 
ences may be due in part to the influence of tip bluntness. It is true 
that the Reynolds nimbers based on tip diameter were higher in experi- 
ments of reference 2 than in the present one. The former were mostly in 
the range 65 < < 230, whereas all the data from the Low Density 

Wind Tunnel correspond to < 9. However, it seems unlikely that 

tip bluntness could account for much of the difference. The experiments 
on flat plates (refs. 10 and ll) indicate that below about Retj_ = 80 

the plate can be considered as sharp, aad the effect of tip bluntness 
has been shown to be much less pronounced for cones than for flat 
plates (ref. 12) . 
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Comparisons Between Theory and Experiment 

Three methods for calculating self- induced pressiires have been em- 
ployed in this report. The first, a new method devised by the present 
authors and called the TC method, is presented and discussed in appen- 
dix A. The second, a method proposed by Probstein, is reviewed and com- 
mented on in appendix B, as is the third, which herein is called the 
TCoo method. Since the differences between the three methods are dis- 
cussed in the appendixes, only the comparisons between theory and ex- 
periment are noted here. 

It can be seen in figures 4 and 6 that the new method for calc\ilat- 
ing self- induced pressures generally overestimates by about 10 to 20 per- 
cent the data obtained in the low-density wind tunnel. In contrast, the 
TCj„ method and the Probstein method using two terms in the series both 
give values greater by a factor of about 2 than the experimental re- 
sults. The better agreement obtained with the present method is not 
surprising; of the three methods it is the only one which accounts in 
even an approximate way for the true interaction effect, wherein the 
^ changes in the external flow due to the presence of the bo\mdary layer 

o feed back into the layer and alter its rate of growth. 

It may be remarked that the present theory has also been compared 
with recent measurments of self- induced pressures on flat plates (ref. 
lO). The agreement was found to be quite good. The term in the expres- 
sion for involving the derivative dM 2 /dx (see appendix A) is not 

negligible in the flat-plate case. When it is included, the calculated 
values of P 2 /P 1 increased by nearly 20 percent in the region near 

the leading edge. The effect of this term was of the order of 5 percent 
at most for the present case, and it was neglected. 

In figure 5, the resilLts of the present theory are represented by a 
single straight line, and the good agreement between the present theory 
and the low-density -wind- tunnel data is clearly evident. _(Actually, the 
individual curves of figs. 4 and 6 when plotted against deviated 

about ±2 to 4 percent from this mean line, but were very nearly linear. 
The deviations do not seem to follow any particular trend and are most 
probably the result of inacc\oracies which were introduced in the graphi- 
cal parts of the analysis.) 

One rather interesting featirre of the data is that they appear to 
be quite linear in X^. Also, the present method gives theoretical 
ctirves which are much more nearly linear than either of the two other 
methods. It may be that the interaction effect results in a boundary- 
layer growth which by coincidence gives induced pressures which are al- 

. most linear in X , at least for X^, < 4. 

c-^ u 
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Hole-Size Effect 


It has been shown by Talbot (ref. l)^ Rayle (ref. ll), and others 
that the apparent pressure sensed by a static -pres sure orifice increases 
with the diameter of the orifice. The phenomenon is due to mixing be- 
tween the stream passing over the surface and the fluid confined within 
the orifice and pressure tubulation. The momentum transferred by the 
mixing sets up currents in the fluid within the orifice which give rise 
to the increase in pressure. 


Ideally^ a static-pressure orifice should be as small in diameter 
as possible^ both to minimize this hole-size effect and to provide a 
truly localized pressure measurement. However^ in rarefied gas flow if 
a pressure orifice is made small enough one encounters another effect 
known as thermal transpiration which can also result in errors in pres- 
sure measurement. The thermal -transpiration effect occurs^ for example, 
when an orifice whose diameter is small compared with the mean free path 
in the gas separates two regions of gas at different temperatures (ref. 
14) . It is easily shown that in this case the pressure ratio is given 

■by 


P2 



( 1 ) 


The static -pressirre orifices used in the present experiments were 
0.010 inch in diameter. For an orifice of this size, the pressure in- 
crement due to momentum mixing is completely negligible. However, there 
is the possibility that thermal-transpiration effects may be important, 
since the boundary layer is a region of strong temperature gradient, and 
many of the molecules which enter the orifice from the gas stream come 
from regions in the boundary layer which are at temperatures different 
from the gas within the orifice. 


A rough estimate of the magnitude of the thermal-transpiration ef- 
fect can be made in the following way. Assume the boundary -layer char- 
acteristics to be given with sufficient accuracy by Howarth*s analysis 
(ref. 15). For an insulated cone, with a = 1 and [- 1/^2 = 



2 

1 + 0.2 M 2 


(2) 

(3) 


Now/ also assume that the velocity distribution in the boundary layer is 
linear in y. Then the temperature at a distance from the wall 

is given by 
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to 

O) 


_ 1 + 0-2 m |(1 - A ^/ S ^) 

^aw 1 + 0.2 m| 

where is the mean free path of the gas at the wall. 

As a specific example, take the following conditions: M 2 = 5.5; 

Re 2 = 10 , 000/ inch; P 2 = 170 microns of mercury; x = 1 inch; and = 

270° K. For these conditions » 0.011 inch and 6 * 0.099 inch. 

Then from equation (l), p^/p^ = 0.994, where p^ is the press\ire in 
the bo\andary layer and p^ is the pressirre within the orifice, which 
is presumably different from py^ because of the thermal-transpiration 
effect. It is seen that for this particular case the error is about I /2 
percent. For the worst conditions encoimtered in the tests, the error 
is about 2 percent. Actually, this analysis greatly overestimates the 
effect, since equation (l) is true only for A^^/d» 1. For the experi- 
ments A^/d “ 1, and in this range the pressure increment is less than 
10 percent of what it is in free molecule flow. One may conclude, 
therefore, that the pressures measured in the experiments were true sta- 
tic pressures, essentially uninfluenced by either momentum mixing or 
thermal-transpiration hole-size effects. 



CONCLUSIONS 

The results of the investigation of viscous self- induced pressures 
on 3° -semi vertex- angle cones may be summarized as follows: 

1. The new data obtained for self-induced pressures on 3° cones are 
in good agreement with prevxous data obtained on 5 cones. All these 
data are correlated reasonably well by the hypersonic-viscous-interaction 

parameter , where M^ and Re^j^ are Mach number and 

Reynolds number based on ideal Taylor-Maccoll flow conditions and C is 
the Chapman -Rube sin factor. 

2. The tangent-cone viscous-interaction model predicts self- induced 
pressures which are only 10 to 20 percent higher than the measured values, 
if local values of the flow parameters at the edge of the boundary layer 
are used in calculating the boundary- layer growth. 


University of California, 

Berkeley, Calif., J^lne 24, 1957. 
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APPENDIX A 


NUMERICAL-GRAPHICAL METHOD FOR CALCULATING SELF- INDUCED 
PRESSURES ON CONE SURFACE 

The self-induced pressure generated on the surface of a cone in hy- 
personic flow because of the outward displacement of the external flow 
streamlines by boimdary- layer growth has been studied by several authors 
(refs. 16 and 17). The assumptions usually made are: (l) A region of 

inviscid flow exists between the outer edge of the boundary layer and the 
shock wave, and (2) the flow parameters, such as Mach number or pressure, 
at the edge of the boundary layer (subscript 2 in fig. l) can be ob- 
tained to satisfy approximations by the tangent-cone (TC) method. The 
tangent -cone method consists of relating the local flow parameters on a 
body to the undisturbed flow ahead of the shock wave through conical-flow 
theory (e.g., Taylor -Mac coll values), but using the local body inclina- 
tion as the effective cone angle. In the case of self- induced pressures 
on a cone as shown in figure 1, the local effective cone angle is taken 

to be 0^^ of the cone angle 0 ^ and the angle 0g = tan"^ 

Two computations are involved in utilizing the TC method for evalua- 
tion of self- induced pressures. First, the inviscid- flow values must be 
obtained, for given effective cone angle 62 ^ either from the exact solu- 
tion (Kopal*s tables) or by one of several approximate methods (refs. 18 
and 19) which are available. Secondly, the boundary-layer displacement 
thickness must be evaluated as a function of position along the cone sur- 
face, and here again several methods of varying accuracy and complexity 
are available. It will be noticed that the two computations are not in- 
dependent, since the boundary-layer growth determines the inviscid- flow 
values at the outer edge of the boundary layer, but at the same time the 
rate of growth of the boundary layer is detemnined by these inviscid- flow 
values. An accurate application of the TC method must include this in- 
teraction effect. 

For the computations of the inviscid- flow values, it was decided 
to use the exact Taylor -Maccoll results as computed by Kopal. There are 
several supersonic and hypersonic flow approxjjnations in analytic form 
which are accurate over different ranges of the similarity parameter 
K 2 = ^±^ 2 ’ single one is sufficiently accurate over the entire 

range of K 2 encountered in the tests reported here. 

For the computation of the growth of boundary- layer displacement 
thickness d5*/dx, it would be desirable to have a method which takes 
into account the external-flow pressure gradient, the transverse- 
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curvatiire effect^ an accurate viscosity- tempera tiore relation^ the appro- 
priate value for the Prandtl number^ and the actual measiured wall tem- 
perature. The formula for 5* which is used here^ and taken from ref- 
erence 20, fulfills all the foregoing requirements with the exception of 
the consideration of the external-flow pressure gradient and transverse- 
curvature effect. The approximate expression of reference 20 employs 
the Chapman-Rubesin factor C evaluated in a way which gives close 
agreement with Crocco*s exact calculation based on the Sutherland law. 
The effect of isothermal wall temperature different from the adiabatic 
value and the effect of Prandtl number are also Included, but the re- 
sults are valid only for zero pressLire gradient. To acco\mt in some 
measure for the effect on the boundary layer of the variation in the 
external-flow quantities along its outer edge, local values of the 
external-flow parameters are used to calculate Mg and ^^^xg* This 

method does not take into account changes in the velocity, temperature, 
and density distributions within the boundary layer caused by the 
external-flow gradients. 


Method of Computation 

The computation method employed is a semigraphical one. 

Inviscid flow; TC approximation . - For a given free-stream condi- 
tion (Mj_, pj^, Tj^, and Re^^), plots of Mg, P 2 /pq^ ^ 2 /"^!^ ^^2 

constructed as functions of equivalent effective cone angle 6g. In- 
terpolation in Kopal*s tables and the adiabatic law are used to obtain 
these values. Viscosities are tsLken directly from a plot of the semi- 
empirical fomula of reference 21. Curves of the form shown in sketches 
(a) and (b) are obtained. 
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Boundary- layer displacement thickness , - The approximate result of 
reference 20, to vhich the Mangier transformation correction l/ -\/3 
(ref. 22) has been applied, is 


X 


dS-^ 

dx 


where 







1 + 





(Al) 


^ = 1 + <r^/2 (r - 1) „2 

To ^ f 2 2 


(A2) 


and the Chapman -Rubes in factor C is defined by 

^ “ tTt^ 

The viscosity i-i' is evaluated at the intermediate temperature T' 
given by 

T'/Te = (VTe) - 0.468 - (T^^/T^)] - 0.273 a IL^ 


(A3) 




(A4) 


However, instead of using the Sutherland law, as recommended in refer- 
ence 20, the values of reference 21 are used for since the 

Sutherland law is less accurate at the low free-stream temperatiires en- 
co\mtered in the present tests. The actual wall temperature Tv is 
determined by experiment. 


For a given free-stream condition (M-j^, p^^, etc.) a set of values of 
02 is chosen and for each of these values a range of values of 0g is 
computed using local free-stream conditions given by plots similar to 
sketches (a) and (b). The angle 0 q is presented as a function of x 
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for each selected value of 62 and may be represented graphically as 
shown in sketch (c). 




(d) 


The value of x appropriate to a particular value of ©2 that 
for which ©2 “ “ ®c^ the cone angle. Thus, the values displayed in 

sketch (c) are replotted with x as a parameter in the form shown in 
sketch {d) . The particular values of © 2 (x) for which ^2 “ " ®c 

are then used in connection with sketch (b) to obtain P 2 /P 1 ^ func- 
tion of X, where P 2 /P 1 cone free-stream pressirre ratio, in- 

cluding the self- induced pressure. The ideal Taylor-I^ccoll pressure 
ratio Pq/pi^ corresponding to ©2 = can of course be obtained di- 
rectly from sketch (b), or in the case of the 3° cones, from Van Dyke's 
second-order theory. 


Discussion 

Accuracy of TC method for calculating inviscid flow . - The TC method 
has been utilized in its most accurate form, namely, with exact values 
taken from the Kopal tables. This was necessary because neither the 
hypersonic approximation of reference 18 nor the Karman-Moore slender- 
body result (ref. 23) is sufficiently accurate for K 2 of the order of 
unity, which is the range encountered in the present tests. However, it 
must be remembered that the TC method itself provides only an approxima- 
tion to the inviscid-flow field. The acc\iracy of the method has been ex- 
amined (refs. 17 and 24) by comparing results obtained by it for pressure 
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distributions on pointed ogives with exact values obtained from the method 
of characteristics. The TC method yields surface pressures which are 
slightly higher than the exact values^ the difference depending on the 
distance from the vertex of the ogive. At the vertex the two methods^ 
of course^ give identical resiilts; farther back the deviation may be of 
the order of a few percent. It seems reasonable to expect the same to 
hold true for the present application^ though a precise correspondence 
cannot be made because the effective shape formed from the cone plus the 
boundary layer is a blunt body. 


It is noted that reference 2 gives an incorrect interpretation to 


the correction factor of reference 17 


A = 


exact 

(dp/ds),pQ ^ 


where 


dp/ds 


represents the initial press\ire gradient at the vertex of a pointed con- 
vex ogive. The factor A is greater than unity, which means the pres- 
sure falls more rapidly back from the vertex according to the exact cal- 
CTolation than it does according to the TC method. Thus the TC method 
overestimates the pressure. In reference 2 the factor is written incor- 


rectly as 


(dp/d0) g^t 

(dp/ d0)ip^ 


and from this it is inferred that the method 


underestimates the pressiare. 


Induced- pressure flow model . - The basic idea used in the TC calcu- 
lation for self- induced press\ires is that the effect of boundary-layer 
growth on the external inviscid flow can be evaluated approximately by 
increasing the effective local cone angle by the amoimt 

= tan"^(d8*/dx) . Order-of-magnitude arguments concerning the accuracy 
of this approximation are given in references 11, 16, and 17. The error 
involved in replacing the actual streamline inclination in the external 
flow by the slope is estimated to be of the order (&/x)^, where 5 

is the bo\mdary-layer thickness. It also tirrns out that the neglect of 
pressiore gradient across the boundary layer is justified provided (5/x)^ 
is small. 


Boundary- layer calculations . - The boundary-layer calculations used 
in this report accotint in some measiire for the variation in external- flow 
properties along the outer edge of the boundary layer by using local val- 
ues for Reg, I>Z’ forth. However, the expression for d5*/dx 

(eq. (A1)) is only approximate, since the terms involving dMg/dx, 
dTg/dx, and dEeg/dx have been neglected. It turns out that for the 
present calciilations these terms contribute an increment of about 5 per- 
cent at the most to 0g , and their neglect is not serious . 

Two effects which have not been accoimted for in the boundary-layer 
analysis are the effect of transverse curvatiire and the direct effect of 
the self- induced pressure gradient on the density and velocity distribu- 
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to 

I 

Eri 

O 


tions within the boundary layer. Both the transverse curvature and the 
pressure gradient tend to thin the boundaiy layer and thus result in 
smaller values for the induced-pressure increment. The transverse- 
curvature effect has been studied by Probstein and Elliott (ref. 25). 
Probstein concludes (ref. 26) from his analysis ^ which is valid for small 
that transverse curvature does not appreciably alter the boimdary- 

layer displacement thickness. However^ much of the present data were in 
the range 6*/r^ =1-3^ and for these values the transverse -curvature 

effect almost certainly cannot be neglected. 


In reference 27 a correction for transverse curvature is deduced which 
/.^ , where is the displacement thickness^ in- 

® Vl + 2S*/re ^ 

eluding transverse curvature, and 5* is the displacement thickness on the 
cone in the absence of curvature effects (i.e., eq. (Al) of this report). 
This result is based on the ass\JH 5 )tion, suggested by the calculations of 

oo 

reference 24, that the integral / ( 1 - = hi is unaffected by 

transverse curvature. Then the identification 6* = ^ 


d6^ 


_r_ 

r^ 


dr is 


If one 


made^ and the result for cited above follows immediately. 

applies this transverse curvature correction to the present data it is 

seen that for — - 1 (roughly the lower limit of the data) the result 
^c 

cited above predicts a decrease in induced pressure of about 40 percent, 

while for the largest value, — = 3, the decrease is about 60 percent. 

^c 

These decrements in pressiore seem much too large, and it is probable that 

d^ 

the assim^tions used to obtain are not good approximations for val- 

ues of — as large as those of the present tests. 

^c 


Evaluation of resiilts of calculations . - It can be seen from figures 
4 (f) and 6 that the method given in this appendix yields induced-pressiare 
increments which are smaller by factors of about 1.5 to 1.8 than those 
obtained from the Probstein and TC methods. The difference is due solely 
to the interaction process, which has been approximated by us in^ local 
values for Mg, Reg, forth in the calculation of dS /dx. 


An improvement in the accuracy of the present method co\ild probably 
be made by accoiuiting for some of the neglected effects which have been 
mentioned in the foregoing discussion. For example, the calculated in- 
duced pressure would be increased by (l) including the gradient terms 
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(3M2/'ix, and so forth, in the expression for d5*/dx, and (2) correcting 
for the error involved in the approximation 0g = tan“^(d5*/dx) . On the 
other hand, the calculated induced pressures would be decreased by: (l) 

including the transverse -cijrvature effect, (2) correcting for the error 
in the TC approximation for the inviscid-flow field, (3) including the 
pressure gradient across the boimdary layer, and (4) Including the direct 
effect on boundary-layer growth of the pressure gradient along the outer 
edge of the bo\mdary layer. 
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APPENDIX B 


PROBSTEIN'S ANALYSIS FOR SELF- INDUCED-PRESSURE EFFECT 


Probstein (ref. 26) considers a Taylor series expansion of the sur- 
face pressure in the form 



Now^ for < 1, which is the range encountered in the present tests ^ 
the Lees hypersonic approximation is not accurate, so that the values r 
for the derivatives in equation (B2) obtained by Probstein cannot be 
used. Probstein suggests that for < 1 the ratio Pp/p^ be evalu- 
ated from the Kopal tables, but that the derivatives be evaluated from 
the Karman slender-body result (ref. 23). However, this resiilt is also 
not sufficiently accurate. Following Van Dyke's suggestions (ref. 28) 
for combined supersonic -hypersonic similarity, a formula of the form 


^=1 

Pi 


AiTMfe^ 

+ 5 log. 



(B3) 


was used to fit the Kopal values and the Van Dyke second-order- theory 
values for Pc/p^ over the range 0.14 < 0^ - 1 < 0.3, 0 < 0g < 0.13 

radian. The constants foiind were A^^ = 1.52 and Ag = 2.85. Ejqoression 
(B3) was used to calculate the derivatives^in equation (B2) , and expres- 
sion (A1) was used to evaluate 6g and 0g, except that inviscid-flow 

values M , Re„, and T were used instead of the local values M 2 , Rs 2 ^ 
and Tg. 
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For 


(B4) 


As a check on the above values^ a numerical differentiation was per- 
formed on a curve constructed from cross plots of the Kopal entries and 
the 3^ values calculated by second-order theory^ and the values 3.8 and 
38 were obtained for the derivatives. The first of these values may be 
more accurate than that given in equations (B4) . Not much accuracy can 
be claimed for either set^ however. 

The results of the Probstein analysis are shown in figures 4(f) and 
6 . The first-order curves were obtained using only the first term in 
the series (B2)^ second-order c-urves using two terms. It can be seen 
that for small-angle cones (i.e.^ small K^) the convergence of the ser- 
ies is slow. This slow convergence is already evident at = 1; as 
can be seen from examination of the functions presented in reference 26. 

As a check on the Probstein analysis^ two induced-pressure distri- 
butions were evaluated by what is called here the TC^ method. In this 
method the boundary-layer slope 65 is evaluated as above, using 
inviscid-flow values of T^, and Re^ (i-e., the values which obtain 

far downstream in the flow field external to the boundary layer), and is 
added point by point to 0^ to obtain 62 us a function of x. The 

pressure distribution P 2 /Pq is obtained directly from a plot of the 
Kopal entries. 

The TC„ method of course gives the values that the Probstein method 
should converge to. As is evident in figures 4(f) and 6 , the use of two 
terms in the Probstein series provides a fairly good approximation. But 
the method of Probstein is not very useful for < 1 because of the 
difficulty in determining the necessary derivatives with any accuracy. 


It was necessary to include at least two terms in the series. 


example, with M 2 _ = 3.70, 0^ = 5 


Pi r 
L 


^(Pc/Pl) 
^0 




= 3.5 


0 = 0 , 


^ p^(Pe/Pl)1 

P- L J 




= 17 


0 = 0 , 
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TABLE I. - FLOW COITOITIONS OF TESTS 
[stagnation temperature, 300° K. ] 


Mach 

number^ 

Ml 

Static 

pressure^ 

Pl^ 

microns 

Hg 

Reynolds 
number 
per inch^ 
Re 

Pressure 

ratio^ 

Pc/Pi 

(a) 

3.70 

49.8 

880 

1.106 

3.91 

73.0 

1,510 

1.114 

3.97 

85.1 

1,860 

1.117 

4.05 

108.8 

2,540 

1.120 

5.47 

66.5 

4,530 

1.207 

5.73 

113.3 

8,980 

1.223 


^Cone semivertex angle 


e^, 3,03°. 



Figure 1. - Schematic drawing of viscous flow over cone. 
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-DIAM. 



L 

D 

X, 

X2 

^3 

X4 

Xj 

^6 

^7 

A 

MODELS 

9.50 

0.500 

0.75 

1.25 

2.00 

2.75 

3.50 

4.00 

4.50 

B 

MODELS 

10.75 

0.750 

0.75 

1.25 

3.00 

4.50 

6.00 

6.50 

7.00 




& 


Vg O.D. SOFT- 
COPPER TUBING 


Figure 2. - Model specifications. Material, all brass except for l/2-inch-outside- 
dlameter soft-copper tube; four static-pressure orifices 0.010 inch in diameter 
located at 90° Intervals on each model at distances given in table (dimensions 
in Inches ) . 
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(a) Free-stream Mach number, 3.70; free-stream Reynolds 
number per inch, 880. 

Figure 4. - Induced-pressure Increment plotted against x. 
Cone semivertex angle, 3^. 


28 


NACA TN 4327 



X, in. 


(b) Free-stream Mach number, 3.91; free-stream Reynolds 
number per inch, l,5io. 


Figure 4 . 


Continued. 
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(c) Free-stream Mach number, 3.97; free-stream Reynolds 
number per inch, 1,850. 


Figure 4. 


Continued . 
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X, in. 


(d) Free-stream Mach number, 4,05; free-stream Reynolds 
number per Inch, 2,540. 


Figure 4. 


Continued. 
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(e) Free-stream Mach number, 5.47; free-stream Reynolds number per 
inch, 4,530. 
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Figure 4. 


Continued. 


MCA TW 4327 



.6 


O) 

CO 



(f) Free-stream Mach number, 5.73; free-stream Reynolds number per 
inch, 8,980. 


Figure 4. - Concluded. 
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Figure 5. - Induced-pressure Increment plotted against viscous-interaction 
parameter Xc* 
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NACA - Langley Field, Va. 





Figure 6. - Induced-pressure increment plotted against x for 
present data and data of reference 1. 
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